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We discuss the properties of Lorentzian wormholes in dilatonic Einstein-Gauss-Bonnet theory 
in four spacetime dimensions. These wormholes do not need any form of exotic matter for their 
existence. A subset of these wormholes is shown to be linearly stable with respect to radial pertur- 
bations. We perform a comprehensive study of their domain of existence, and derive a generalised 
Smarr relation for these wormholes. We also investigate their geodesies determining all possible 
particle trajectories, and perform a study of the acceleration and tidal forces that a traveler crossing 
M ' the wormhole would feel. 
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I. INTRODUCTION 



Black holes are by far the most celebrated class of solutions derived from Einstein's field equations. Being 
among the first types of solutions to be found almost a century ago, they have undergone an extensive investi- 
gation over the years. Their existence conditions, forms of solutions and set of properties have been studied in 
the context of both traditional General Relativity and generalized gravitational theories admitting either four 
or more dimensions. 

The most well-known example of such a generalized gravitational theory in four dimensions is provided by 
, the low-energy heterotic string effective theory [HQ- In this theory, the scalar curvature term of Einstein's 
theory is only one part of a more complex action functional where higher-curvature gravitational terms as well 
\ as kinetic and interaction terms of a variety of additional fields (axions, fermions and gauge fields) make their 

■ appearance. The dilatonic Einstein-Gauss-Bonnet theory is a minimal version of the aforementioned theory 

■ and contains the scalar curvature term R, a scalar field called the dilaton, and a quadratic curvature term, the 
Gauss-Bonnet (GB) term, given by i?^^ = R^upaR'^'^'"^ —ARf^^R'^'^ +R^. The GB term in four dimensions can 
be expressed as a total derivative term and, normally, makes no contribution to the field equations; however, 
the exponential coupling to the dilaton field, that emerges in the context of the dilatonic Einstein-Gauss-Bonnet 
theory, ensures that the GB term is kept in the theory. 

In the framework of the dilatonic Einstein-Gauss-Bonnet theory, new types of black hole solutions emerged 
that are endowed with non-trivial dilaton hair (for an indicative list of works on this topic, see The 
' presence of the GB term in the theory caused the circumvention of the traditional no-hair theorems as it 
bypassed the conditions for their validity. In reality, the existence of the dilatonic black holes violated only 
the 'letter' of the no-hair theorems, and not the 'essence' of them since the dilaton charge was of a 'secondary' 
type. However, it became evident that in the context of this type of generalized gravitational theories solutions 
with a much richer structure and a modified set of properties, compared to the ones in General Relativity, can 
emerge. 

Another class of gravitational solutions whose properties are strictly set by the General Theory of Relativity 
are wormholes. They were first discovered in 1935 as a feature of Schwarzschild geometry Q and named the 
'Einstein-Rosen bridge' as they connect two different universes. Their importance was realised in full when 
Wheeler [Tol . [Tl| showed that such a bridge, or wormhole, can connect also two distant regions of our own 
universe thus opening the way for fast interstellar travel. However, it was soon demonstrated that this is not 
possible for the following reasons: (i) the Schwarzschild wormhole is hidden inside the event horizon of the 
corresponding black hole, therefore, is not static but evolves with time; as a result, the circumference of its 
'throat' is not constant but opens and closes so quickly that not even a light signal can pass through [l^ [l3|. 
(ii) even if a traveler could somehow pass the throat, she would be bound to exit the wormhole through the 
past horizon of the Schwarzschild geometry; this horizon was shown to be unstable against small perturbations 
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and that it would change to a proper, and thus non-traversable, horizon at the mere approaching of the traveler 

[iia. 

In hj a new class of wormhole solutions was found that possess no horizon and could in principle be 
traversable. However, some form of exotic matter whose energy-momentum tensor had to violate all (null, weak 
and strong) energy conditions was necessary in order to keep the throat of the wormhole open. Several studies 
have considered a phantom field, i.e. a scalar field with a reversed sign in front of its kinetic term, as a candidate 
for such a form of matter |l8l-[23| . 

It was demonstrated in [J] that the GB term leads to an effective energy-momentum tensor that also violates 
the energy conditions. It is in fact this violation that causes the circumvention of the no-hair theorem forbidding 
the existence of regular black holes with non-trivial scalar hair. It was noted recently in [2^, too, that the 
presence of the GB term in the context of a scalar-tensor theory has the property to evade the various no-go 
theorems of General Relativity. Indeed, various wormhole solutions were found in the context of gravitational 
theories with higher curvature terms |24l - l27| . In the presence of the Gauss-Bonnet term in particular, wormhole 
solutions were found in the context of higher-dimensional gravitational theories [28l - [36| . 

In this work, we will investigate the properties of wormhole solutions that arise in the context of the four- 
dimensional dilatonic Einstein-Gauss-Bonnet theory, first reported in [37l] . The presence of the higher-curvature 
GB term, that follows naturally from the compactification of the 10-dimensional heterotic superstring theory 
down to four dimensions, suffices to support these types of solutions without the need for phantom scalar fields 
or other forms of exotic matter. 

The outline of our paper is as follows: In section II, we present the theoretical context of our model and 
discuss the asymptotic forms of the sought-for wormhole solutions at the regions of radial infinity and the 
regular throat. Based on the latter, we derive the embedding diagram and study the violation of the energy 
conditions. We also present a Smarr relation for the wormhole solutions. In section III, we present the results of 
our numerical analysis that reveal the existence of wormhole solutions in the dilatonic Einstein-Gauss-Bonnet 
theory, and discuss their properties. We demonstrate the stability with respect to radial perturbations of a 
subset of these solutions in section IV. In section V we discuss the junction conditions. The geodesies in these 
wormhole spacetimes are presented in section VI. We calculate the magnitude of the acceleration and tidal 
forces that a traveler traversing the wormhole would feel in section VII, and conclude in section VIII. 



II. EINSTEIN-GAUSS-BONNET-DILATON THEORY 



A. Action 

We consider the following effective action 0, [s^ motivated by the low-energy heterotic string theory [H 
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where <j) is the dilaton field with coupling constant 7, a is a positive numerical coefficient given in terms of the 
Regge slope parameter, and R^^ = R^vpaR^^'"^ — ^R^.uR'^" + R^ is the GB correction. 
The dilaton and Einstein equations are given by 
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with 



(2) 



(3) 



H^y — 2 ^RR^y — 2R^pR^ — 2R^pjya-R'^^ + Rfipo-xRjf^^l ~ 2^^^^gb ' 



(4) 
(5) 



3 



B. Ansatz and equations 



Throughout this paper we consider only static, spherically-symmetric solutions of the above set of equations. 
Thus we can write the spacetime line-element in the form [J] 



ds^ = g^^ (Ixf^dx" = -e^'-^'^dt^ + e^'-^'^dr^ + {d9^ + sin^ Odip^) 



(6) 



As was demonstrated in the dilatonic-Einstein-Gauss-Bonnet (EGBd) theory admits black hole solutions 
whose gravitational background has the line-element of Eq. It was also shown that further classes of solutions 
emerge in the context of the same theory. One of them, in particular, possesses no curvature singularity and 
no proper horizon (with gu being regular over the whole radial regime). However, the grr metric component as 
well as the dilaton field showed some pathological behavior at a finite radius r — vq, as seen from the expansion 
near ro [1] 



= Xi{r-ro) + 
r'(r) = ■ 



t>o + 4>Wr - ro + ■ 



(7) 
(8) 
(9) 



The absence of any singular behaviour of the curvature invariants at ro signifies that the aforementioned patho- 
logical behaviour is merely due to the particular choice of the coordinate system. 

In [s^l we have argued that this class of asymptotically flat solutions can be brought to a regular form by 
employing the coordinate transformation r^ = Vq + P. Then, the metric becomes 



ds^ ^ -e^-i'^dt^ + f{l)df + {f + rl) {de^ + i 



(10) 



The above form is regular and describes a wormhole solution, where ro is the radius of the throat. Indeed, in 
terms of the new coordinate I, the expansion at Z = assumes the form 



fil) = h + hi + 



= 00 + 01^ + 



(11) 
(12) 
(13) 



where fi, Vi and are constant coefficients, and shows no pathology. Both metric functions and the dilaton 
field remain finite in this asymptotic regime. Thus these solutions possess no horizon. In addition, all curvature 
invariant quantities - including the GB term - turn out to be finite at ^ = 0, a result that demonstrates the 
absence of any singularity. 

Substitution of the metric Eq. (fTO|) into the dilaton equation ([2]) and Einstein equations ([3]) yields a coupled 
system of ordinary differential equations (ODE's) for the metric functions and the dilaton field 
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Here r'^ = P + Tq and the prime denotes the derivative with respect to I. Equations and (|16l) fohow 

from the tt, II and 99 components of the Einstein equations, respectively, whereas the last equation (|T7)) follows 
from the dilaton equation. For the numerical computation we 'diagonalize' Eqs. (ITU) . ([T5|) and ([TT]) with respect 
to /', v' and 4>" ■ The remaining equation, Eq. (|16p. involving also second derivatives of u and is satisfied if 
the other three equations are fulfilled. Thus a system of ODE's must be solved that consists of two first-order 
equations for the metric functions / and u and a second-order equation for the dilaton field. 



C. Expansions 

The expansion near the throat, Eqs. (ITT|) - ([T^ . once substituted into the set of equations leads to a number 
of recursive constraints that determine the higher-order coefficients in terms of the lower ones. Looking at the 
lowest order, we observe that /o, vq and 4>q are free parameters. Also, the set of parameters of the theory 
includes the radius of the throat rg and the value of a. The value of the constant 7 will be later fixed to 1 for 
simplicity. 

However, not all of the above parameters are actually independent. To start with, we observe that the field 
equations remain invariant under the simultaneous changes 

c^^<i> + c^., {r,l) ^ {r,l) e~^'f"/\ (18) 

In addition, the following transformation 

In k 

a^ka, (j)^(f)-\ . (19) 

7 

is also a symmetry of the equations. In the light of the above, we conclude that only one parameter out of the 
set (a,ro,<5!io) is independent. We may therefore choose to have a zero asymptotic value of the dilaton field at 
infinity which entails fixing the value of the dilaton field at the throat ^o- The remaining two parameters can 
be combined to give a dimensionless parameter a/rQ that will be used throughout our analysis. Finally, among 
the two parameters associated with the metric functions, {i^q, /q), again only the latter is independent - it is 
only the derivatives of the metric function that appear in the field equations of motion, therefore, we may use 
this freedom to fix the value of in order to ensure asymptotic flatness at radial infinity. Thus, our class of 
wormholes is a two-parameter family of solutions that, in the context of our analysis, have been chosen to be 
/o and a/rg. 

When the expansion of the metric functions and dilaton field near the throat pT|l - (|13p are substituted into 
the field equations, we obtain a set of constraints on the higher-order coefficients. The constraint on the value 
of the first derivative of the dilaton field at the throat is particularly interesting and takes the form 

-bl = /o(fo-^^ _ . (20) 
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As the left-hand-side of the above equation is positive-definite, the same must hold for the right- hand- side. We 
may easily see that the expression inside the square brackets in the denominator remains positive and has no 
roots if 
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'0 ^ 

This inequality is automatically satisfied for the set of wormhole solutions presented in the next section, and 
gives a lower limit on the size of the throat of the wormhole tq . Then, the positivity of the right-hand-side of 
Eq. (PO)) demands that fo > 1. The solutions satisfying /o = 1 comprise a boundary in the phase space of the 
wormhole solutions and their physical significance will be discussed in the next section. 

As already discussed above, at the asymptotic regime of radial infinity, i.e. as ^ — > c», we demand asymptotic 
flatness for the two metric functions and a vanishing value of the dilaton field. Then, the asymptotic expansion 
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at infinity takes the form 

M 

+ (22) 
D 

^ ^ -- + ... . 

In the above, M and D are identified with the mass and dilaton charge of the wormhole, respectively. We 
remind the reader that in the case of the black hole solutions 4], the parameters M and D were related, and 
that rendered the dilatonic hair as "secondary" . However, in the case of the wormhole solutions these parameters 
are not related; this result, together with the fact that the number of independent parameters near the throat 
is also two, confirms the classification of this group of solutions as a two-parameter class of solutions. 



D. Wormhole geometry 



A general property of a wormhole is the existence of a throat, i. e. a surface of minimal area (or minimal 
radius for spherically symmetric spacetimes). Indeed, this property is implied by the form of the line element 
([T0| above, with /(O) and i'(O) finite. To cast this condition in a coordinate independent way, we define the 
proper distance from the throat in the following way 



If we impose the condition for a minimal radius R = -^/P + Tq at Z = 0, this translates to 

(PR 



dR 
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It is easily seen that the first condition is indeed satisfied. For the second condition we find 
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This gives a coordinate independent meaning to the parameter /o. 



/o 



R 



fR 



(23) 



(24) 



(25) 



(26) 



A particularly useful concept with which we may examine the geometry of a given manifold is the construction 
of the corresponding embedding diagram. In the present case, we consider the isometric embedding of a plane 
passing through the wormhole. Due to the spherical symmetry of the solutions, we may simplify the analysis 
and choose 6 = Tr/2. Then, we set 



f{l)dl^ + {P + rl)dip^ ^ dz"^ + dif + ri^dip^ 



(27) 



where {z,r],ip} are a set of cylindrical coordinates in the three-dimensional Euclidean space R^. Regarding 
and 7] as functions of Z, we find 



From the last equation it follows that 

z{l) = ± 



= fil) 



dl' 



(28) 



(29) 
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Hence {ri(l), z{l)} is a parametric representation of a slice of the embedded 9 = 7r/2-plane for a fixed value of 
(p. We observe that the curvature radius of the curve {r]{l), z{l)} at Z = is given by Rq = ro/o- Thus fo has a 
geometric meaning: fo = Ro/rg is the ratio of the curvature radius to the radius of the throat. 



E. Energy conditions 

As was already mentioned in the Introduction, the existence of the wormhole solution relies on the violation 
of the null energy condition. The null energy condition holds if 

T^^nf^n" > (30) 

for any null vector field n''. In the particular case of spherically symmetric solutions, and if we employ the 
Einstein equations, this condition can be expressed as 

-G°o + G\>0, and - GjJ + > . (31) 

If one or both of the above conditions do not hold in some region of spacetinie, then the null energy condition 
is violated. By making use of the expansion of the fields near the throat (fTT|) - (IT3|) . we find that, close to rp, 

[-G"o + G\]^^^ = --^<0, (32) 

provided e^-^^o) ^ 0, i.e. in the absence of a horizon. Thus for the wormhole solutions there is always a region 
close to the throat where the null energy condition is violated. 

On the other hand, from the asymptotic expansion of the solutions at infinity, we find 

-G° + Gl ^ ^i + 0{l-') , (33) 

-G'^o + Gt ^ 20aMD^ + O{l-'^) . (34) 

We observe that, for solutions with a positive dilaton charge D (and a positive mass M) the null energy condition 
is satisfied in the asymptotic region. However, if the dilaton charge is negative, the null energy condition is 
violated also in that asymptotic region. Note that D > for the black hole solutions, where it was given 
by a positive-definite combination of {a, M,(f>oo) [!]• However, D > does not necessarily hold for all of the 
wormhole solutions. 



F. Smarr relation 



To derive the Smarr-like mass formula for the wormhole solutions we start with the definition of the Komar 
mass 

M = Mth + R^^.e^n^dV = Mth " ^ / RoV^d^x , (35) 

where is the timelike Killing vector field, E is a spacelike hypersurface, n'^ is a normal vector on S and dV 
is the natural volume element on E. Here Mth denotes the contribution of the throat, 

Afth = l^th^ : (36) 
where Ath is the area of the throat, and k is the surface gravity at the throat 

Ac2 = -l/2(V^e.)(V^r), K^^iy'iO). (37) 

v/o 
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To obtain the mass formula we express Rq in terms of the effective stress energy tensor 
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where we have added a zero in the form of the dilaton equation. MultipUcation by \/—g yields 
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Substitution of the metric yields a total derivative for the right-hand-side of the above equation. Thus, after 
integration, we find 
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Now substitution into Eq. (|35|) gives the Smarr-like formula 



M = 
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Here h is the induced spatial metric on the throat, R is the scalar curvature of h, and the integral is evaluated 
at Z = 0. Defining the normal vector on the surface I — the mass formula becomes 
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with 



Dth = ^ J Vhe'">n^df,(j) (l + 2a-f'^e-'<'f'R^d?x 
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(44) 



According to the above, the Smarr-like mass formula for wormholes is obtained by replacing the horizonprop- 
erties by the corresponding throat properties in the known EGBd mass formula for black holes [1, |^. In 
addition, an extra contribution appears that may be interpreted as a modified throat dilaton charge, while the 
GB modification is of the same type as the GB modification of the area (or entropy in the case of black holes) . 



III. NUMERICAL WORMHOLE SOLUTIONS 



For the numerical calculations we use the line-element (jlOp . since the functions are well behaved at rp. For 
the representation of the numerical results we employ the metric as well. 

To solve the ODEs numerically, we introduce the compactified coordinate x — + I), thus mapping the 
semi-infinite range of I to the finite range of x, i.e. < x < 1. We cover the parameter range 0.002 < a < 0.128 
and 1.0001 < fo < 20.0, keeping ro ~ 1 fixed. Also we set 7 = 1. 



A. Metric and dilaton functions 

Let us start the discussion of the solutions by recalling the boundary conditions for the system of equations, 
consisting of two first order and one second order ODEs. At the throat I — regularity requires 
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FIG. 1: (a) The metric function /(/), (b) the metric function v(l\ (c) the dilaton function and (d) the scaled GB 
term a7?QB versus I for wormhole solutions with the parameter values /o — 1.1, 1.5, 2.0 and a/rg — 0.02, 0.05. 



This boundary condition has to be supplemented by /(O) = /o or 0(0) — (f>o in order to obtain a specific 
solution. We note that the asymptotic condition / 1 for Z oo is always satisfied. This can be seen from 
the asymptotic form of Eq. (14), /' + /(/ — l)/l = 0, which has the general solution f = + const). This 
leaves the asymptotic boundary conditions at radial infinity 

lim = , lim (/) = . (46) 

Wormhole solutions can be found for every value of a/rg below a/rg w 0.13. This upper bound on a/r^ 
translates into a lower bound on the radius of the throat tq (for a given a). Thus the radius of the throat vq 
can be arbitrarily large. 

In Figs. la,b and c, we show the metric functions f{l) and v{l) and the dilaton function 0(/), respectively, 
for an indicative set of wormhole solutions. We also exhibit the scaled GB term aR^^ in Fig. Id. Close to the 
throat the functions show a distinct dependence on both parameters, /o and a/r^. However, for intermediate 
values of Z, the functions f{l), and likewise the functions ^{l), corresponding to the same value of a/rg, approach 
each other and form clusters. For larger values of I, solutions obtained for different values of a/r^ also merge 
together. The same behaviour is observed for the GB term whereas for the dilaton function (p{l) the same 
tendency also exists but becomes pronounced at slightly larger values of I. 

In Fig. [2] we visualize the geometry of the wormhole solutions. In particular, we present as a typical example 
in Fig. [2^ the isometric embedding of the solution with a/rg = 0.02 and /g = 1.1 [33|- Here also the curvature 
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FIG. 2: (a) Isometric embedding of the wormhole solution for a/rg — 0.02 and fo = 1.1 (taken from [s^l), (b) z versus 
rj for a/r^ = 0.02 and fo = 1.1, 2.0, 20.0. 

radius at the throat, Rq = tq/o, is shown. Also for a/r^ — 0.02, the curves z versus i] near the throat are shown 
in Fig. [2Jd for several values of fo. 

We may also define the profile function b{r) via the equation 

e-Hr) ^ 1 _ M . (47) 
r 

We note that, at the throat, the profile function b/ro goes to one and, thus, g^"^ vanishes. For small values of 
fo, the value of the metric function v at the throat keeps decreasing as /o — >■ 1. Therefore, in the limit /o — >■ 1, 
the metric function gtt tends to zero and a horizon emerges, thus recovering the class of dilatonic black hole 
solutions [1]. 

B. Domain of existence 

Let us now explore the domain of existence of these wormhole solutions. To that end we construct a set of 
families of wormhole solutions, where for each family the value of the parameter a/rg is fixed, while the second 
parameter fo varies within a maximal range of 1 < /o < oo. The values of the parameter a/rp cover the range 
< a/rQ < 0.13. For larger values of a/r^ no wormhole solutions are found. 

In Fig. [3] we present the domain of existence of the wormhole solutions. Fig. [5^ shows the scaled area of the 
throat A/16ttM'^ versus the scaled dilaton charge D/M. The domain of existence is mapped by the families 
of solutions obtained for a representative set of fixed values of a/rQ. We observe that the domain of existence 
of the wormhole solutions is bounded by three curves indicated by asterisks, crosses and dots. The boundary 
indicated by asterisks coincides with the EGBd black hole curve [J| and corresponds to the limit fo I. The 
limit /o — > oo on the other hand is indicated by crosses. At the third boundary, marked by dots, solutions are 
encountered that are characterized by a curvature singularity. 

Before discussing the three limiting cases in more detail, let us consider further the boundary conditions at 
the throat. We illustrate the metric function fo at the throat for the same set of families of wormhole solutions 
in Fig. [SJd. For fixed values of a/r^ below the limiting value for EGBd black holes, a/?'o|]-,h = 0.0507 fo 
covers the full range 1 < fo < oo. Beyond this critical value, however, the minimal value of fo increases with 
increasing a/r^. We note that, for a certain intermediate range of a/rp, fo is not monotonic. 

In Fig.[3l: we show the dilaton field (j)o at the throat versus the scaled dilaton charge D/M for the same set of 
solutions. Clearly, the domain of 0o is bounded. Again we observe that in the limit /o ^ 1 the black hole values 



10 




-9 ^ ' ' ' ' 1 

-1 I • • ' • ' • • • • ' 1 0.2 0.4 0.6 0.8 1 

-0.1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1 « 



(c) (d) 

FIG. 3: (a) Scaled area A/IGtyM^ of the throat versus the scaled dilaton charge D/AI for several values of a/r-Q. The 
boundaries represent EGBd black holes (asterisks), limiting solutions with /o — >■ oo (crosses) and solutions with curvature 
singularities (dots), (b) Metric function /o at the throat (determining the curvature radius Ro — ro/o) versus the scaled 
dilaton charge D/M; (c) dilaton function <^(0) at the throat versus the scaled dilaton charge D/AI; (d) redshift function 
i/{0) at the throat versus l//o- 

are obtained. Inspection of the redshift function i^q at the throat, exhibited in Fig. [5Ji, reveals that —500(^0) 
tends to zero as /o — > 1. Thus a horizon emerges in this limit. 

1. Black hole limit 

To study the black hole limit in more detail let us consider a sequence of wormhole solutions approaching 
the black hole solution. We demonstrate this limiting behaviour for solutions with a = 0.02 in Fig. |4l Here, 
we exhibit a sequence of solutions with values of /o tending to one. Clearly, the profile function b{r) (Fig. 2^) 
and the dilaton function (/){t) (Fig. [4)d) tend fast to the limiting black hole solutions. For /o = 1.0001 the 
wormhole functions and their black hole counterparts are already very close. The redshift functions v are not 
distinguishable for this set of solutions except very close to the throat, where diverges in the limit. 



11 





(a) 



(b) 



FIG. 4: (a) Profile function 6/ro and (b) dilaton function (p versus r/ro for fo 
togetfier with tlie limiting black hole functions. 



1.1, 1.01, 1.001, 1.0001 and a/rg = 0.02 



2. Large fo limit 



20 I.,.,-- ' -. j. 'i — ^ — — — — ^ — ^ 1.8 
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(a) (b) 

FIG. 5: (a) Metric function f{l) and (b) scaled GB term aR^^ versus I for fo = 10.0, 15.0, 20.0 and Of/rg = 0.02, 0.1. 

Next we consider the limit of large fo, indicated by crosses in Figs. [3l As /o — >■ oo, the mass M and the 
dilaton charge D assume finite values. The same holds for the redshift function i/q and the dilaton field 0o at 
the throat, although the derivative of the dilaton field (pQ with respect to the coordinate I diverges like ^/fo, as 
seen from the boundary condition Eq. (jlSj). 

As an example we exhibit the function /(Z) in Fig. [5^ for two values of a/rQ and increasing values of /q. 
We observe that for a given value of a/rQ, the functions f{l) deviate from each other only for small values of I 
close to the throat, but coincide for larger I. Moreover, the region where the functions coincide increases with 
increasing values of /q. Thus the solutions approach a limiting solution for /g — >■ oo, which depends on a/r^. 

This limiting behaviour is even more pronounced in the scaled GB term uRq^, a curvature invariant which 
is exhibited in Fig. [5)d for the same set of parameters. Extrapolating the value of ai?QB at the throat for large 
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/o indicates that ai?QB remains finite for /□ — >■ oo. 



3. Singularity limit 



s 



(a) 



100000 



10000 



(c) 




(b) 



(d) 




FIG. 6: (a) Scaled GB term a7?GB 

versus I, (b) metric function /(/), (c) derivative f'{l), and (d) derivative of the dilaton 
function (f)' {I) for several values of /o and a/rp = 0.1 demonstrating the emergence of a curvature singularity. 



Finally we turn to the third boundary curve, indicated by dots in Figs. [31 This boundary emerges when 
branches of solutions with fixed a/rg terminate at singular configurations. We demonstrate that these singular 
configurations possess a curvature singularity at a critical value of Z, Zcrit, in Fig. [6^. Here the scaled GB term 
ai?Q3 is shown versus I for a/r^ = 0.1 for a sequence of solutions approaching the singular configuration. We 
observe a sharp peak in the vicinity of Zcrit, which increases in size and diverges as the singular configuration is 
approached. 

To elucidate this emergence of a curvature singularity at some finite critical value oil, we show in Fig. the 
metric function /(Z), in Fig. [Bt its derivative f'{l), and in Fig. [6ji the derivative of the dilaton function (j)'{l) 
for a fixed value of a/r^ = 0.1 and increasing values of /o approaching the singular configuration. We note 
that all functions are continuous, but as the singular configuration is approached, their derivatives develop a 
discontinuity at some point lent outside the throat. In fact the derivative of the function / has a pole at Zcrit- 
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(a) (b) 
FIG. 7: Demonstration of the violation of the nuU energy condition for ajr^ — 0.02 and several values of /q. 



C. Energy conditions 

We now turn to the energy conditions, Eqs. ([5T|) . We demonstrate the violation of the null energy condition 
in Figs. [T] Here we show 'normalized' quantities in order to emphasize the change of sign. The normalization 
factor is 



(48) 



This normalization is responsible for the steep rise observed in the figures for /o = 1.1, which occurs when all 
quantities are close to zero. 

The null energy condition is violated close to the throat for all wormhole solutions, as seen explicitly from 
Eq. (|32|) and demonstrated for the set of solutions of Fig. [TK- However, for large values of /o and at the same 
time small values of ajr^ the dilaton charge D can become negative. In this case the null energy condition is 
violated also in the asymptotic region, as seen from Eq. (l34t . In fact we observe that for the larger values of 
/o in Fig. [TId the combination — Tq -f Tq is negative everywhere. Thus in this case the null energy condition is 
violated in the whole spacetime. 



D. Smarr relation 



Let us finally turn to the Smarr- like mass relation, Eq. (j43|) . This mass relation is a perfect check for the 
numerical accuracy and thus for the quality of the solutions. We demonstrate the Smarr relation in Figs. [5] 
by comparing the mass M obtained from the asymptotic fall-off of the metric functions with the mass Msmarr 
obtained by evaluating the various terms in the Smarr relation Eq. (j43p . 

The Smarr-like formula is well satisfied for all sets of solutions. In Fig.[8K, we keep 7 = 1 as in the main body 
of the paper and exhibit the masses M and Msmarr versus I//0 for ajr^ = 0.1. The relative error in this case 
is below 3 • 10"^ (for ajrl = 0.05 it is below 10-^ and for a/rg = 0.02 it is below 5 • 10"^). 

In Fig.lHlD we have fixed = 0.02, /o = 1.4 and varied 7 to address the 7-dependence of the mass formula 
for a set of solutions. Also for these solutions the relative error is small and, in particular, it is below 4 • 10~^. 

To gain a better understanding of the quantities contributing in the Smarr formula, we now consider the 
various throat properties. First, we exhibit in Fig. |9K the scaled throat quantity S'th/47rM^ versus the scaled 
dilaton charge D/M for the full domain of existence of the wormhole solutions. S'th/47rM^ resembles in its 
structure the entropy of black holes, since it contains the same correction term to the throat area ^th as in the 
black hole case to the horizon area. In the black hole limit, ^th assumes the meaning of the black hole entropy. 
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The scaled entropy-analogue throat quantity S'th/47rM^ is always greater than one. Only in the Schwarzschild 
black hole limit, where the dilaton vanishes, it assumes precisely the value one. A comparison with Fig. [3^ 
shows, that the quantity StYi/lQirM'^ is always larger than the scaled throat area A/1QttM^ . This was observed 
before for black holes ^ Q. The color coding of Fig. [3^ is related to the stability of the solution, discussed in 
the next section. 

The scaled surface gravity ktq at the throat of the wormhole solutions is exhibited in Fig. [5]d. It is bounded, 
with its lower boundary given by the black hole values, which remain close to the Schwarzschild value of 1/2. 
The product of the quantities S'th and k represents the first term of the Smarr formula, up to numerical factors. 
We exhibit this first term, divided by the mass, in Fig. |9]:. Since the first term is always greater than one 
(except for the Schwarzschild case), the remaining pieces of the mass formula must contribute negatively to 
cancel this excess. 

The second term in the mass formula, is negative, except for the small region with large /o and small 

a/rp. Its contribution to the scaled mass formula can be read off the horizontal axis of Fig. [Q}:. We exhibit 
the third term of the scaled mass formula in Fig. IHJi. It contains the scaled throat dilaton charge, I?th, which 
is modified by the GB term, analogous to the modification of the area by the GB term. As expected, this 
contribution is negative for the wormhole solutions. It vanishes for the black holes, which represent the upper 
limit for Z3th- 



IV. STABILITY 



Our starting point for the study of the stability is the line element for spherically symmetric solutions 

ds^ = -e^^dt^ + fdl^ + {l^ + rl) {dd^ + sin^ ddip^) . (49) 

Here we consider only the pulsation modes. Thus we allow the metric and dilaton functions to depend on the 
radial coordinate I and the time coordinate t, 

D = u{i,t), f^fii.t). 4>^ki.t)- (50) 

For the study of the stability behaviour of our solutions, we also need the time-dependent Einstein and 
dilaton equations. These were presented in in the context of the stability analysis of the dilatonic black 
hole solutions, and thus we refrain from repeating them here (43 |. Next, we decompose the metric and dilaton 
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FIG. 9: (a) Scaled entropy-analogue throat quantity S'th/47rM^ versus the scaled dilaton charge D/M for several values 
of a/rQ. The boundaries represent EGBd black holes (asterisks), limiting solutions with /o — )■ oo (crosses) and solutions 
with curvature singularities (dots). The shaded areas indicate linear stability (lilac or lower), instability (red or upper), 
undecided yet (white) w.r.t. radial perturbations (Sect. IV). (b) Scaled surface gravity tiro at the throat versus the scaled 
dilaton charge D/M; (c) scaled first term of the mass formula; (d) scaled third term of the mass formula. 



functions into the unperturbed functions and the perturbations 

i){l,t) = v{l) + €5v{l)e"'^ , (51) 
f{l,t) = f{l) + e6f{l)e'''' , (52) 
^il,t) = + eS(j}{l)e"'* , (53) 

where we assume a harmonic time dependence of the perturbations and e is considered as small. Now we 
substitute the perturbed functions into the Einstein and dilaton equations and linearize in e. This yields a 
system of linear ODEs for the functions 5v{l), 5f{l) and 54){l), where the coefficients depend on the unperturbed 
functions and their derivatives. We use the unperturbed equations to eliminate v' , v" , f and 0". 

The tl part of the Einstein equations yields Sf in terms of Scf) and Scf)' . Thus df and 6f' can be eliminated 
from the rest of the ODEs. The dilaton equation and the 96 part of the Einstein equations can be diagonalized 
with respect to Sv" and 6(j)" . Finally Sv' can be eliminated from the dilaton equation by adding the rr part of 
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the Einstein equations with a suitable factor. Thus we end up with a single second order equation for Scj), 

5^" + qM' + (qo + qa(T^)6<j) ^ , (54) 

where the coefficients qi, qo and q^ depend on the unperturbed solution. 

The coefficient q^ tends to one asymptotically and is bounded at ^ = 0. The coefficients qi and qo tend to 
zero asymptotically, however, they diverge at Z = as l/l. Thus, to obtain solutions which are regular at ^ = 
suitable boundary conditions are needed. Since the perturbations are assumed to be normalizable, 6(1) has to 
vanish asymptotically. In order to obtain a unique solution we also fix Scf) at I ~ 0. Thus we are left with three 
boundary conditions for a second order ODE, which can be satisfied only for certain values of the eigenvalue 
a^. 

We can avoid the singularity of qi and qo a,t I = 0, by employing the transformation Scf) = F{l)il){l), where 
F{1) satisfies F' / F = -qi{l)/2. This yields 

tP" + Qo^ + a^q,^ = , (55) 

where the coefficient Qo = ^Qif^ ~ 9i/4 + qo is bounded at Z = 0. We note that F{1) diverges like l/l at I — 0. 
Thus for acceptable solutions we have to impose the boundary condition tp = at I = 0. In addition ijj has 
to vanish asymptotically to ensure normalized solutions. However, these boundary conditions do not lead to 
unique solutions, since the ODE (|55p is homogeneous. Therefore, we supplement Eq. (|55p with the auxiliary 
ODE TV' = tfi'^ and impose the boundary conditions iV(0) = and N 1 asymptotically. These conditions 
give exactly the normalization ip'^dl — 1. Again, solutions exist only for certain values of the eigenvalue cr^. 

The ODE (|55p can be written as a Schrodinger equation by introducing the new coordinate y via dy/dl = ,Jq^, 
and eliminating the first derivative dip/dy as above. This yields 



fx 

dy 



KffX + a^X = , (56) 



where 



1 



9<t/ 4 



1 



2 



and X is related to tp. The new potential Vee is bounded for all y. 

Note that the derivation of Eq. ([55)1 holds only if the function q^r is positive for all I. However, ii q^ < 
on some interval, one cannot rely on the rule that the smallest eigenvalue corresponds to the eigenfunction 
without knots. Indeed we found in some cases two different eigenfunctions ipi, ■02 without knots which satisfy 
the orthogonality condition J qa-'ipi'ip2dl = 0. 

In order to determine the change of stability we consider families of wormholc solutions with fixed ct/rQ. The 
ODE ((54l) is then solved together with the normalization constraint for varying values of /q. In Fig. [TOT a) we 
show the eigenvalue ct^ versus /o for several values of a/r^. We observe that solutions for negative values of 
exist for large values of /o, corresponding to an instability of the wormhole solutions. For these solutions 
the function 64>, respectively ^, decays exponentially. Thus these solutions correspond to bound states of the 
equivalent Schrodinger equation. As /q decreases, the eigenvalue increases and tends to zero for some critical 
value fo — fo"\ which depends on a/rp. No solutions were found for fo < f^"'\ Solutions of Eq. ^ with 
positive would oscillate in the asymptotic region, and are not normalizable. 

The interpretation is that the instability mode vanishes at some critical value /q . Hence we conclude that 

the wormhole solutions are stable for fo < fo ■ However, this holds only if > for all I. Otherwise one 
cannot be sure that is the smallest possible eigenvalue. Therefore, eigenfunctions with negative eigenvalues 
cannot be excluded even if ct^ = 0. Consequently, if is not positive for all I the question of stability cannot 
be decided by considering the standard equivalent Schrodinger eigenvalue problem. 

Numerically, we found that q„ is positive for all wormhole solutions with a/rg < 0.05; but wormhole solutions 
where q„ is negative on some interval exist if a/rg > 0.05. In Fig. llOf b) we show the stability region in the 
domain of existence [s^l • Also shown is the line where the eigenvalue cr^ changes sign. The instability region 
(red or upper) is characterized by negative , whereas the stability region (lilac or lower) is characterized by 
the absence of negative eigenvalues and positive qa- In the remaining (white) region the question of stability is 
not yet decided. 
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(a) (b) 

FIG. 10: (a) Eigenvalue cr^ versus /o for a/r^ = 0.02. (b) Domain of existence: scaled area A/IQ-kM^ of the throat 
versus the scaled dilaton charge D/M for several values of a/rg. The shaded areas indicate linear stability (lilac or 
lower), instability (red or upper), undecided yet (white) w.r.t. radial perturbations (taken from [37|). 



V. JUNCTION CONDITIONS 



Up to this point, we have discussed the behaviour of the metric functions and of the dilaton field in only half 
of the wormhole spacetime, i.e. the part with I > 0. Our solutions should naturally be extended to the second 
asymptotically flat part of the manifold {I — > —oo). If this is performed by demanding that the derivatives of 
the metric and dilaton functions are continuous, we observe a singular behaviour corresponding to curvature 
singularities. This is demonstrated in Fig. [TT]for the wormhole solution with parameters a/rg = 0.02, /o — 1.1 
and 7=1. 
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FIG. 11: Non-symmetric and symmetric extension of (a) the metric function / and (b) the dilaton function for 
a/ro — 0.02 and /o = 1.1 versus I. 



However, wormhole solutions without curvature singularities can be constructed when we extend the wormhole 
solutions to the second asymptotically flat part of the manifold in a symmetric way. In this case jumps appear 
in the derivatives of the metric and dilaton functions at I — 0. The jumps can be attributed to the presence of 
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matter located at the throat of the wormhole. The corresponding junction conditions are of the form 



0:76 



Sdil 



(57) 



where sj^ is the stress-energy tensor of the matter at the throat, and Sdii the corresponding source term of the 
dilaton field. The Ihs of the junction conditions can be derived in a standard way by integrating the Einstein 
and dilaton equations across the boundary ^ 0, i. e. {G"^ — T^) = ^ lim^^^o J^l i^l^ ~ -^m) '^^^^ yields 



T 



J a 



, 



(G'o-T^) 



, a7e 



-700 



'0 > 



(58) 
(59) 
(60) 
(61) 
(62) 



where the subscript indicates evaluation at Z = 0. 

Next we assume that the matter at the throat takes the form of a perfect fluid with energy density p and 
pressure p and a dilaton charge pdii- We also introduce the action 



(Ai -I- Ao2ae-'"^i?) V-hd^x 



(63) 



at the throat, where hab denotes the (2-|-l)-dimensional induced metric on the throat, R the corresponding 
Ricci scalar, and Ai, A2 are constants. Inserting the metric this brings the non-trivial junction conditions 
Eqs. (I5ll)-(IS21) to the form 



8a7e' 



P — Ao 2 '^i ) 



8a7e" 



-700 



To 



= A, 



—700 



Pdil 

2 



(64) 
(65) 
(66) 



Using these equations p, p and pdii can be expressed in terms of the metric and dilaton functions and the 
constants Aq, Ai. In Figs. [T2]we give an example for Ao = Ai = 1 and 7 = 1. We note that the stable wormhole 
solutions possess positive energy density p. 

Let us also consider the special case p = (i.e. dust) and choose the dilaton charge density pdii at the throat 
to be twice the dilaton charge density of the wormhole in the 'bulk', pdii = '^4'o/Vfo- This yields 



Ai — Ao — 



2^0 



(67) 



8a7e' 



-700 



To 



Interestingly, Ao — Ai and p are positive for all wormhole solutions, as shown in Fig 1131 

Next we consider perturbations of the wormhole solutions in the form f{l,t) = /(/) + (5/(Z)e*'^*, etc. For 
symmetric wormholes the derivatives of the perturbations also have a jump at E. Introducing perturbations of 



19 




-0.8 -0.6 -0.4 -0.2 



(a) 




-0.8 -0.6 -0.4 -0.2 



0.2 0.4 0.6 0.8 



1 



-e- 
I 

a 



>.Q=1.0, ^,=1.0 




-0.8 -0.6 -0.4 -0.2 



1.7 
1.6 
1.5 
1.4 
1.3 
1.2 
1.1 





'^0 



(b) 



),n=1.0, X, = -l.O 



-0.8 




0.002 
0.01 
-I 0.02 
0.03 
0.04 
0.05 
0.052 
0.055 
0.06 
0.07 
0.08 
0.09 
-I 0.10 
0.11 
0.12 
0.125 
0.128 
stab 



0.2 0.4 0.6 0.! 



0.002 
0.01 
0.02 
0.03 
0.04 
0.05 
0.052 
0.055 
0.06 
0.07 
0.08 
0.09 
0.10 
0.11 
0.12 
0.125 
0.128 
stab 



(d) 



FIG. 12: (a) Energy density p , (b) pressure p, (c) dilaton charge density pdii and (d) interaction potential versus 4>o- 



the energy density Spe^"* , pressure Spe''"^ and dilaton charge (5pdiie"^* we find 



5p 
6p 
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(69) 
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(71) 



Thus these equations give the perturbations of the energy density, the pressure and the dilaton charge on 
the throat. Thus for any perturbation in the 'bulk', we can find the corresponding perturbation on the throat. 
Consequently, the stability of the solutions is not affected by the matter on the throat. 
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FIG. 13: (a) Energy density p and (b) A = Ao = Ai versus 00 for dust. 

VI. GEODESICS 



The study of the orbits of test particles and hght is essential to fully understand the properties of a spacetime. 
The motion of test particles in a EGBd wormhole spacetime with metric 5^^, Eq. ()10|) . and dilaton field (j) is 
governed by the Lagrangian 



(72) 



where the dot denotes the derivative with respect to the afhne parameter r, and /3 is a constant (/3 — 1/2 for 
heterotic string theory). The conjugate momenta = are found to be 



= -e -'Pfe'-^i , PI = e'^P^fl , pg = e'^^^r^ + 1^)6 , p^ = g-^^^r^ + f) sin^ 



First integrals of motion are given by 

Pt = const. 



-E , Pip = const. = L 



(73) 



(74) 



We refer to E as the energy of the test particle and to L as its angular momentum. 

The affine parameter can be chosen such that 2C = k, with k — — \ for time-like geodesies and k — Q for null 
geodesies. Since we are considering spherically symmetric spacetimes we may choose = 7r/2, which implies 
= 0. Employing Eqs. (|74p the Lagrangian then reduces to 



2£ = e'P^e-^" 



Let us first consider time-like geodesies. For k = — 1 we obtain 



K . 



/ 



where we introduced the effective potential 



V3,(l,L) = e'-'(^ke 



[E'-K%{l,L) 



.-2/30 



(75) 



(76) 



(77) 
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The effective potential Kff is a suitable quantity to discuss qualitatively the trajectories of test particles. We 
note that the effective potential approaches the value of one asymptotically. Thus it follows from Eq. (|76p that 
unbound trajectories are only possible if _E > 1. The turning points of trajectories for a given energy E and 
angular momentum L are determined by the condition 

E = V,sir,,L). (78) 

We note that, in contrast to black hole spacetimes, the wormhole spacetimes do not possess an event horizon. 
Consequently, particles cannot disappear behind an event horizon. Instead they can travel through the throat 
of the wormhole from one asymptotically flat part of the manifold to the other. 
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FIG. 14: Effective potential Vch versus logio(''/''o) for several values of the angular momentum L for the wormhole 
solution with a/r-Q = 0.05 and /o — 1.1. 

To discuss the various types of trajectories in these wormhole spacetimes we consider as an example the 
wormhole solution with parameters a = 0.05 and /o = 1.1. We exhibit the corresponding effective potential 
Vef[(r, L) in Fig.[T3]for several values of the angular momentum L. In particular, we first discuss the dependence 
of the shape of the effective potential on the angular momentum L of the test particle, and next consider the 
different types of trajectories with their dependence on the energy E of the test particle. 

In the case L — Q the effective potential is a monotonically increasing function of r, assuming its minimum at 
the throat r — tq and tending to one asymptotically. When the angular momentum is increased, the effective 
potential remains monotonic until it develops a saddle point for some critical value L^j-n. For L > icrit the 
effective potential is no longer monotonic, since a local maximum Vma.x{L) and a local minimum VminiL) occur. 
If L is larger than a certain value, Li (say), the maximum T4nax(i) exceeds the asymptotic value of the effective 
potential, and thus it becomes the global maximum. We note that the essential features of the effective potential 
Kff (r, L) are the same for all the wormhole solutions considered. 

Concerning the types of orbits of test particles we note that there are two kinds of trajectories. For trajectories 
of the first kind the particles remain on a single asymptotically flat manifold, whereas for trajectories of the 
second kind the particles travel from one asymptotically flat part of the manifold to the other, passing through 
the throat. 

Trajectories of the first kind exist only if the effective potential is non-nionotonic, i.e. for L > La-it- Bound 
orbits exist for Vniin(-^) < E < 1, e.g. E2 in Fig. [T31 whereas unbound orbits exist for 1 < i? < Kiiax(i), e.g. E3 
in Fig. 1141 For unbound orbits we need L > Li, thus VmaxiL) must be the absolute maximum. 

Trajectories of the second kind exist for all values of E above the minimum of the effective potential Vcs{rQ, L). 
The particles move on bound orbits if VcB{ro,L) < E < max{l, Vmax(i)}, e.g. Ei, E2, E^ in Fig. [141 Thus 
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for E = Ves{ro,L) the particles move on circles at the throat. Radial trajectories are obtained for L = 0. 
In this case the particles oscillate about the throat. In the limit Eg — y/Vcs{rQ, 0) they possess the minimal 
possible energy and are at rest at the throat. On the other hand, ii E > max{l, Vmax(i)}, e.g. E4 in Fig. [TH 
the particles travel along unbound orbits, starting at infinity on one asymptotically fiat part of the manifold, 
passing through the throat and reaching infinity on the other asymptotically flat part of the manifold. 

In Fig. [15] we exhibit the minimal possible energy Eq = y^Vcs{ro, 0) versus D/M for several values of a/r 



Note that Eq approaches zero in the black hole limit, 
indicates the stability change of the wormhole solutions. 



The black line on the left, highlighted in the inset, 
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FIG. 15: Minimal energy Eo = ^/Vcsifo, 0) of test particles versus the scaled dilaton charge D/M for several values of 
a/rQ. The black line on the left denoted by cr'^ = indicates the stability change of the wormhole solutions. 



In order to calculate the trajectories we consider 



dl I r?i + 



from which we find 



(p{l) ^ ipo±L 



lo '0 



/ 



r2 + P \ E^ - V: 



-dl' 



(79) 



where I is restricted to intervals where E^ — V^g > and tpo is an integration constant. The trajectories are 
then displayed in the xy-plane, where 



= ^rl+Pcosiip{l)) , y{l) = y^r^ + P sin(^(Z)) . 



(80) 



In Fig. [16] we show examples of bound orbits of test particles with angular momentum L — 2 and several 
values of the energy E for the wormhole solution with a/rp = 0.05 and fo = 1.1. For an energy of E = Ei = 0.9 
(Fig. 116b ) there exists only a bound orbit of the second kind. This is highlighted in the figure by using different 
colors (line styles) for the two asymptotically flat manifolds, when projecting into the ccy-plane. In contrast, for 
E — E2 — 0.98 (Fig. \Wb) there exists in addition a bound orbit of the first kind. Examples for unbound orbits 
are shown in Fig. [T7k for energies E = E^ = 1.02 and E — E4 = 1.05. 

Let us now consider the null geodesies in these wormhole spacetimes. Setting k = in Eq. (j75p we find 



E' 



E' 



(81) 
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a) 



- 



b) 






FIG. 16: Bound orbits of massive test particles with angular momentum L — 2 and energy (a) E — 0.9 and (b) E = 0.98 
in the wormhole spacetime with a/rg = 0.05 and fo — 1.1. The colours red (solid) and blue (dotted) indicate motion on 
the first and the second asymptotically flat part of the manifold, respectively. The throat of the wormhole is shown by 
the black circle. The right panel of a) shows the orbit on the isometric embedding of the wormhole. 



We define an effective potential Vcs{l, L/E) by 

Kla,Vi^)-(|)'^^(|)^eyO. (82) 

We note that for the discussion of the null geodesies it is more convenient to consider the "normalized" effective 
potential Wcfi(0 = , since it is independent of L/E. 

As an example, we show the "normalized" effective potential v^s for the parameter values a/r^ = 0.05 and 
/o = 1.1 in Fig. [T71 The "normalized" effective potential possesses a local minimum timin at the throat, a 
maximum timax at some distance from the throat and tends to zero asymptotically. 

As for the timelike geodesies there are two kinds of trajectories: either the massless test particle remains on 
one of the asymptotically flat parts of the manifold, or it passes through the throat from one asymptotically 
flat part to the other. The first kind of (unbound) trajectories exists only if the ratio E/L is smaller than the 
maximum of the "normaUzed" effective potential. The second kind of trajectories includes unbound geodesies 
if E/L > Wmax and bound geodesies if Vmin < E/L < Wmax- Circular orbits exist for E/L = Vi^in = Vcsil = 0) 
(and unstable ones for E/L — Umax)- 
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FIG. 17: (a) Same as Fig.[T6]for unbound orbits of massive test particles and energies E — 1.02, 1.05. (b) "Normalized" 
effective potential Vcft for massless test particles versus logj^Q(r/ro) for the wormhole solution with a/rg = 0.05 and 
/o = 1.1. 



VII. ACCELERATION AND TIDAL FORCES 



Following the formalism of [17[ , we will now calculate the magnitude of the acceleration and tidal forces that a 
traveler traversing the wormhole would feel. For this purpose, it is particularly convenient to make two changes 
of reference frames: first, starting from the standard reference frame with basis vectors (et, e/, eg, e,^), in which 
the line-element takes the form of Eq. ([TU| , we change to an orthonormal reference frame with 



1 



1 



1 



er = e e 



t, 



-j=ei, eg = -60, 



r sm ( 



(83) 



in terms of which the metric tensor assumes the form: = • = Vap- Alternatively, we may write the 
set of equations in (|83|) as = e^, with 



e-" 

l/Vf 

1/r 
l/(rsin6l) 



The above allows us to write the transformation law of the components of the Riemann tensor as 



net 



^13 ^7 ^ '^P'y ' 



(84) 



(85) 



with L^^ = (L^) ^. Next, we introduce the orthonormal reference frame of the traveler which is related to the 
previous one by a Lorentz transformation 



6(5 , 



(86) 



with 7 = [1 — (u/c)^)] and v = ^{\/Jdl/e'^dt) the radial velocity of the traveler at radius I as measured by 
a static observer there. As before, we may write = A-^ e/i, with 



A^ 



7 Tl{v/c) 

j{v/c) =F7 

10 

1 



(87) 



25 



so that 



^ = A" A-" A-" R'^--- 



We note that • ^ —1, and thus can be naturally considered as the traveler's normalized vector of 
four-velocity u. For the magnitude of u to be fixed in the rest frame of the traveler, the four-acceleration 
a = du/dr should be orthogonal to the four-velocity, i.e. a • u = 0, and therefore a = (0,a') - if we further 
assume that the traveler moves radially, then a = aej. For a traveler moving along their world-line in which 
the tangent vector is = u, the acceleration is given by the formula a — VuU or a'^ = wf^u^c^. In order to 
compute the magnitude of the acceleration that the traveler feels, we work in the following way: we calculate 
the at component in the (t, l,9,Lp) coordinate frame first according to the formula 



c c ^/ J 



where we have assumed that the four-velocity is a function of the radial variable I and rewritten its expression 
as 

V 1 

u = -le'^Bt T 7 75 = it*et + e; . (90) 

However, it also holds that 

V _ V 

at ^ a - et ^ (aej) ■ et = a-f- e "{et ■ et) = -07- e" . (91) 

Combining the above two results for at, we find that the magnitude of the acceleration is 

2 



a 



e-"^ [e-'^y . (92) 



For the wormhole to be traversable, the above quantity must remain always finite and, if possible, take a small 
value. 

In , it was also demanded that the gravitational acceleration should be small at the location of the stations 
where the trip starts and ends. Demanding that the time needed to complete the whole trip should not be too 
large, the stations should be fairly close to the throat. Both of these conditions are satisfied in our case, since 
an asymptotically flat regime is reached fairly quickly. 

The above also means that the traveler, starting from the stations, does not need to travel with a relativistic 
velocity in order to approach the wormhole in a reasonable time. Then, if ?; ^ c, 7 ~ 1, and the acceleration 
reduces to 

H=c^^. (93) 

In Fig. [181 we depict the dimensionless acceleration a = |a|/(c^/ro) at the throat as a function of D/M. 
We observe that a ranges roughly between 10 and 100 for a set of wormhole solutions with parameter a/r^ 
ranging between 0.128 and 0.01. In the context of superstring theory, a ~ ^p, therefore, for the above solutions 
To ~ lO^p and the magnitude of the acceleration turns out to be (10^^ — 10^^) g®, where is the acceleration 
of gravity at the surface of the earth. Since there is no upper limit for rg in our solutions, one may be tempted to 
increase the size of the throat of the wormhole so that |a| is of the order of g^; this would demand a throat radius 
of the order of at least (10 — 100) light-years. The aforementioned result could perhaps have been anticipated 
since our analysis is performed within the context of superstring theory, a theory whose fundamental scale is 
tied to the Planck scale. 

Next we turn to the tidal acceleration that a traveler feels between two parts of her body as the wormhole is 
crossed. Let w be the vector separation between these two parts, which in the reference frame of the traveler 
is purely spatial, w = (0, w*). Then, the tidal acceleration is given by the expression [l^ 
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FIG. 18: The dimensionless acceleration a — \a\/{c^ /ro) that a traveler would feel traversing the wormhole as a function 
of D/M for a variety of values of a/rQ. 



where k takes on only spatial values and where we have used that m° = (5?. But since the metric is diagonal 
and the Riemann tensor is antisymmetric in the first two indices, the superscript p. should also take only spatial 
values. Then, the non-vanishing components of the tidal acceleration are: 



lift 



(95) 



The transformation laws (|84p and (j87p allow us to compute the components of the Riemann tensor in the 
orthonormal frame of the moving observer in terms of the ones in the {t, I, 9, ip) coordinates. We then obtain 



Aa' 



Similarly, we find: 

Aa^ = Aa^ 



' HI 



2 e 2 d9 

-c w I 7 -ft 



tet 
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2 6 2 

-C W J 



r2 + P 



tet 



R' 



wi 



(96) 



(97) 



As in the case of the acceleration, we should demand that the magnitude of the above components should be 
small, i.e. of the order of the acceleration of gravity at the surface of the earth. Using the fact that |w| ~ 2 m, 
the above constraints may be written as 



/V 
2/ 

If 



- v'^ - v" 



P 



< 



< 



9® 



1 



c2(2m) 



9® 



(10i"cm)2 ' 
1 



c2(2m) (10i"cm)2 



(98) 



(99) 



The second inequality involves the velocity with which the traveler moves and thus may be considered as a 
constraint on this quantity. The first inequality restricts again the profile of the metric functions - a similar 
analysis to the one above leads to the same results regarding the magnitude of the tidal forces and the necessary 
size of the throat in order to bring these down to a reasonable value. 
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VIII. CONCLUSIONS 

The existence of traversable wormholes in the context of General Relativity relies on the presence of some form 
of exotic matter. However, in the framework of a string-inspired generalized theory of gravity, the situation may 
be completely different. Here we have investigated wormhole solutions in EGBd theory, which corresponds to a 
simplified action that is motivated by the low-energy heterotic string theory. Indeed, as we have demonstrated, 
EGBd theory allows for stable, traversable wormhole solutions, without the need of introducing any form of 
exotic matter. The violation of the energy conditions, that is essential for the existence of the wormhole solutions 
[43} , is realized via the presence of an effective energy-momentum tensor generated by the quadratic-in-curvature 
Gauss-Bonnet term. 

We have determined the domain of existence of these wormhole solutions and shown that it is bounded by 
three sets of limiting solutions. The first boundary consists of the EGBd black hole solutions of Q. The second 
boundary is approached asymptotically, when the curvature radius at the throat of the wormhole diverges (the 
/o — 00 limit). Finally, at the third boundary, solutions with a curvature singularity at a finite distance from 
the throat are encountered. 

We have investigated the properties of these EGBd wormholes and derived a Smarr-like mass relation for 
them. In this, the horizon properties in the black hole case are replaced by the corresponding throat properties 
of the wormholes, thus the area and surface gravity here refer to the ones at the throat. Moreover, as is well- 
known for black holes in EGB theories, their entropy does not correspond merely to their horizon area but it 
receives a GB correction term. Similarly, we find that the mass formula for the wormhole solutions includes 
an analogous GB correction term; in addition, another term, that vanishes in the black hole case, appears that 
represents the GB corrected dilaton charge at the throat. We have demonstrated that the Smarr relation is 
satisfied very well by the numerical solutions. 

We have also investigated the stability of the solutions. We have shown that a subset of our wormhole 
solutions, the one that lies close to the border with the linearly stable dilatonic black holes in the domain of 
existence, is also linearly stable with respect to radial perturbations. While we have also shown that another 
subset is unstable, we have concluded that the study of the standard equivalent Schrodinger equation cannot 
determine the stability for the full domain of existence We hope to resume the question of the existence of an 
alternative method for the study of the stability of all of our wormhole solutions at some future work. 

When the wormhole solutions are extended to the second asymptotically flat region, this extension must be 
made in a symmetric way, since otherwise a singularity is encountered. As a consequence of this symmetric 
extension, the derivatives of the metric and dilaton functions become discontinuous at the throat. This discon- 
tinuity demands the introduction of some matter distribution at the throat. We have shown that this may be 
realized by the introduction of a perfect fluid at the throat whose energy density is positive for the subset of 
stable wormhole solutions. 

Next, we have studied and classified the geodesies of massive and massless test particles in these wormhole 
spacetimes. Depending on the respective effective potential, there are two general kinds of trajectories for these 
particles. The particles may remain on bound or escape orbits within a single asymptotically flat part of the 
spacetime, or they may travel from one asymptotically flat part to the other on escape orbits, and travel back 
and forth on bound orbits. 

In addition, we have calculated the acceleration and tidal forces which travelers traversing the wormhole 
would feel. We find that their magnitude may be small for fairly large values of the size of the throat of the 
wormhole. According to our findings, the radius of the wormhole throat is bounded from below only, therefore, 
the wormholes can be indeed arbitrarily large. Astrophysical consequences will be addressed in a forthcoming 
paper as well as the existence of stationary rotating wormhole solutions in the EGBd theory. 



28 



Acknowledgments 

We gratefully acknowledge discussions with Eugen Radu. B.K. acknowledges support by the DFG. 



D. J. Gross, J. H. Sloan, Nucl. Phys. B291 (1987) 41. 

R. R. Metsaev, A. A. Tseytlin, Nucl. Phys. B293 (1987) 385. 

D. G. Boulware and S. Deser, Phys. Lett. B 175, 409 (1986). 

P. Kanti, N . E. Mavromatos, J. Rizos, K. Tamvakis, E. Winstanley, Phys. Rev. D 54 (1996) 5049 
•arXiv:hep-th/9511071 . 

Z. -K. Guo, N. Ohta and T. Torn, Prog. Theor. Phys. 120, 581 (2008) |arXiv:0806."248T] [gr-qc]] : Prog. Theor. Phys. 
121, 253 (2009) arXiv;0811.3068 [gr-qc]]. 

N. Ohta an d T. T orn, Prog. Theor. Phys. 121, 959 (2009) arXiv:0902.4072' [hep-th]]; Prog. Theor. Phys. 122, 
1477 (2009) [arXi^: 0908.3918 [hep-th]]; Prog. Theor. Phys. 124, 207 (2010) arXiv:1004. 2779] [hep-th]]. 
K. -i. Maeda, N. Ohta and Y. Sasagawa, Phys. Rev. D 80, 104032 (2009) ,arXiv:0908.4151 [hep-th]]: Phys. Rev. D 
83, 044051 (2011) arXiv: 1012.0568 [hep-th]]. 

B. Kleihaus, J. Kunz, E. Radu, Phys. Rev. Lett. 106, 151104 (2011) [arXiv:1101.2868 ' [gr-qc]]. 

A. Einstein, N. Rosen, Phys. Rev. 48 (1935) 73. 
J. A. Wheeler, Annals Phys. 2, 604-614 (1957). 

J. A. Wheeler, Geometrodynamics (Academic, New York, 1962). 

M. D. Kruskal, Phys. Rev. 119, 1743-1745 (1960). 

R. W. Fuller, J. A. Wheeler, Phys. Rev. 128, 919-929 (1962). 

L H. Redmount, Prog. Theor. Phys. 73 (1985) 1401. 

D. M. Eardley, Phys. Rev. Lett. 33, 442-444 (1974). 

R. M. Wald, S. Ramaswamy, Phys. Rev. D21, 2736-2741 (1980). 

M. S. Morris, K. S. Thorne, Am. J. Phys. 56, 395-412 (1988). 

H. G. Ellis, J. Math. Phys. 14, 104-118 (1973). 

H. G. Ellis, Gen. Rel. Grav. 10, 105-123 (1979). 

K. A. Bronnikov, Acta Phys. Polon. B4, 251-266 (1973). 

T. Kodama, Phys. Rev. D18, 3529-3534 (1978). ^ 

C. Armendariz-Picon, Phys. Rev. D65, 104010 (2002) gr-qc/0201027 . 

K. A. Bronnikov, E. Elizalde, Phys. Rev. D 81, 044032 (2010) a rXiv:0910.3929l [hep-th]]. 

D. Hochberg, Phys. Lett. B251, 349-354 (1990). 

H. Fukutaka, K. Tanaka, K. Ghoroku, Phys. Lett. B222, 191-194 (1989). 
K. Ghoroku, T. Soma, Phys. Rev. D46, 1507-1516 (1992). 

N. Furey, A. DeBenedictis, Class. Quant. Grav. 22, 313-322 (2005) |gr-qc/0410088] . 

B. Bhawal, S. Kar, Phys. Rev. D46, 2464-2468 (1992). 

G. Dotti, J. Oliva, R. Troncoso, Phys. Rev. D75, 024002 (2007) hep-th706 07062l . 

E. Gravanis, S. Willison, Phys. Rev. D75 (2007) 084025 gr-qc/0701152 . 

G. Dotti, J. Oliva, R. Troncoso, Phys. Rev. D76, 064038 (2007) arXiv:0706.11MI [hep-th]]. 

M. G. Richarte, C. Simeone, Phys. Rev. D76 (2007) 087502 arXiv:0710. 2041 [gr-qc]]. 

C. Garraffo, G. Giribet, E. Gravanis and S. Willison, J. Math.^ Phys. 49, 042502 (2008) jarXiv:0711.2992l [gr-qc]]. 

H. Maeda, M. Nozawa, Phys. Rev. D78, 024005 (2008) arXiv:0803.1704 [gr-qc]]. 

C. Garraffo and G. Giribet, Mod. Phys. Lett. A 23, 1801 (2008) arXiY:Q8Q5 3575 [gi-_qc]]. 

C. Simeone, Phys. Rev. D83 (2011) 087503. " 

P. Kanti, B. Kleihaus, J. Kunz, Phys. Rev. Lett. 107 (201 1) 271101 [ar Xiv:1108.3003 [gr-qc]]. 
S. Mignemi, N. R. Stewart, Phys. Rev. D 47 (1993) 5259 [arXiv:hep-th/9212146 . 

C. M. Chen, D. V. Gal'tsov, N. Ohta, D. G. Orlov, Phys. Rev. D81 (2010) 024002 arXiv;0910. 3488 [hep-th]]. 

P. Kanti, N. E. Mavromatos, J. Rizos, K. Tamvakis, E. Winstanley, Phys. Rev. D57 (1998) 6255-6264 

hep-th/9703192 . 

M. H. Daouda, M. E. Rodrigues and M. J. S. Houndjo, Eur. Phys. J. C 71, 1817 (2011) _arXiv:1108.2920l [astro- 
ph.CO]]; Eur. Phys. J. C72 (2012) 1890 |arXiv :1109.0528 [physics.gen-ph]]. 

In the line-element (|6]) was used for the derivation of the time-dependent field equations. The new set of 

equations with respect to the new coordinate / may be derived under the redefinitions F = 2iy, A = ln[r^///^] and 
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[43] According to recent studies [4i[|, wormhole solutions arise also in the context of f{T) gravitational theories, where 
T is the torsion, without the energy conditions being violated. 



